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Abstract 

We extract the strong coupling from the recently updated ALEPH non-strange spectral functions obtained 

from hadronic t decays. We apply a self-consistent analysis method, first tested in the analysis of OPAL data, to 
extract and non-perturbative contributions. The analysis yields - 0.296 + 0.010, using Fixed Order 

Perturbation Theory (FOPT), and = 0.310 + 0.014, using Contour Improved Perturbation Theory (CIPT). The 

weighted average of these results with those previously obtained from OPAL data give = 0.303 + 0.009 and 

- 0.319 + 0.012, which gives, after evolution to the Z boson mass scale, - 0.1165 + 0.0012 and 

= 0.1185 + 0.0015, respectively. We observe that non-perturbative effects limit the accuracy with which 
can be extracted from r decay data. 
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1. Introduction 

The extraction of as from hadronic t decays repre¬ 
sents an important test of the evolution of the strong 
coupling as predicted by the QCD /J-function. At and 
around the r mass, ~ 1.78 GeV, perturbative QCD 
can still be used, but realistic analyses must include the 
contribution from non-perturbative effects. The stan¬ 
dard framework to describe hadronic t decays is to orga¬ 
nize the QCD description in an operator product expan¬ 
sion (OPE) where, apart from the perturbative contribu¬ 
tion and quark-mass corrections, the QCD condensates 
intervene [1]. 
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Observables such as R^, 

_ r[T“ ^ v^hadrons(7)] 

r r - M 

F [t ^ Vre vAr)\ 

can be written as weighted integrals over the experi¬ 
mentally accessible QCD spectral functions. The spec¬ 
tral functions have been determined at LEP by the 
ALEPH [2] and OPAL collaborations [3]. In the spe¬ 
cific case of Rj, the weight function is that determined 
by T-decay kinematics and the integral runs over the to¬ 
tal energy of the hadronic system in the final state, s, 
from zero to m^. Since the OPE description is not valid 
at low-energies the evaluation of the theoretical coun¬ 
terpart is performed exploiting the analytical properties 
of the QCD correlators. One writes then a finite-energy 
sum rule (FESR) where the theoretical counterpart of 
the observable is obtained from an integral along a com¬ 
plex circle of fixed radius |s| = m^. In fact, any analytic 
weight function gives rise to a valid sum rule, and it 
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has become customary to exploit this freedom in order 
to analyse several FESRs simultaneously. This type of 
combined analysis allows for the extraction of as well 
as non-perturbative contributions. 

The computation, in 2008, of the NNNLO term, 
of the perturbative expansion of the QCD cor¬ 
relators [4] triggered several reanalyses of as from t 
decays [5-8]. In the process, it was discovered that 
the correlation matrices of the then publicly available 
spectral functions from the ALEPH collaboration had a 
missing contribution from the unfolding procedure [9]. 
(It was for this reason that we restricted our attention to 
OPAL data in the analyses of Refs. [7, 8].) Recently, 
a new analysis of the ALEPH data became available, 
employing a new unfolding method, which corrects 
for this problem in the correlation matrices [10, 11]. 
The ALEPH spectral functions have smaller errors than 
opal’s and have the potential to constrain the theoret¬ 
ical description better. The new set of ALEPH spectral 
functions motivates the present reanalysis. 

On the theoretical side, two different aspects have re¬ 
ceived attention recently. The first one regards the use of 
the renormalization group in improving the perturbative 
series. There are several prescriptions as to how one 
should set the renormalization scale. The most com¬ 
monly used are Contour Improved Perturbation The¬ 
ory (CIPT) [12] and Lixed Order Perturbation Theory 
(LOPT) [13]. Different prescriptions lead to different 
results with the available terms of the perturbative ex¬ 
pansion, and therefore to different values of as- This 
discrepancy remains one of the largest sources of uncer¬ 
tainty in as extractions from t decays. Strong evidence 
in favour of the LOPT prescription has been given in 
Refs. [14, 15] but the issue is still under debate. Here, 
we chose to perform our analysis using both prescrip¬ 
tions and hence quote two values of a^. 

The second point that has been studied recently is the 
description of non-perturbative effects. Since the work 
of Ref. [5] it is known that the OPE parameters obtained 
in some of the recent analyses are inconsistent. With 
these parameters one cannot account properly for the 
experimental results when sq, the upper limit of the in¬ 
tegration in the EESR, is lowered below m^. A strategy 
that allows for a self-consistent analysis is the inclusion 
of Duality Violation (DV) effects in the theoretical de¬ 
scription. It is well known that in the vicinity of the 
Minkowski axis the OPE alone cannot account for all 
non-perturbative effects. In the past, in the description 
of Rr, this contribution was systematically ignored due 
to the fortuitous double zero of the weight function at 
the Minkowski axis. In the same spirit, combined anal¬ 
yses of several EESRs were restricted to the so-called 


pinched moments, i.e., moments that have a zero at the 
the Minkowski axis. Until recently, DVs were not tack¬ 
led directly and their contribution to final results and 
errors were not systematically assessed. 

Recent progress in modeling the DV contribu¬ 
tion [17-21] has allowed for analyses that include them 
explicitly in the EESRs. In Refs. [7] a new analy¬ 
sis method taking into account DVs explicitly was pre¬ 
sented. This led to a determination of from OPAL 
data in a self-consistent way, together with the OPE con¬ 
tribution and DV parameters [8]. In a recent work, we 
applied the same analysis method to the updated version 
of the ALEPH non-strange spectral functions [22]. 

In the remainder we discuss the main results of 
Ref. [22]. 


2. Analysis framework 


Eor the sake of self-consistency, here we make a brief 
review of the framework of our analysis. The details can 
be found in the original publications Refs. [7, 8, 22]. 

Pits performed in our analysis are based on EESRs of 
the following form [1, 16] 


C(^o) 




where the weight-functions w(s) are polynomials in s, 
and is given by 


/ J d^xe'^’^{Q\T[j^{x)Jl{0)] 10) 

= {q^qy - q^g^y) H- , (3) 


with s - q^ — -0^, and one of the non-strange V or 
A currents or uy^ysd. The superscripts (0) and (1) 
refer to spin. In the sum rule Eq. (2), is the ex¬ 

perimentally accessible spectral function. We construct 
EESRs at several values of sq ^ for ^ given weight 
function. 

The correlators can be decomposed exactly 

into three parts 

(4) 

where “pert” denotes perturbative, “OPE” refers to OPE 
corrections of dimension larger than zero (including 
quark-mass corrections), and “DV” denotes the DV con¬ 
tributions to 
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It is convenient to write the perturbative contribution 
in terms of the physical Adler function [13, 14], that sat¬ 
isfies a homogeneous renormalization group equation. 
When treating the contour integration, one must adopt a 
prescription for the renormalization scale. As discussed 
above, we perform our analysis within CIPT and FOPT 
and quote both results. 

The leading contributions from higher dimensions in 
the OPE can be parametrized with effective coefficients 
Co as 

= (5) 

The dimension-two quark-mass corrections can safely 
be neglected for the non-strange correlators.' There¬ 
fore, in our analysis C 2 = 0.^ The first non-negligible 
contribution is then C 4 , that can be related to the gluon 
condensate. However, the weight functions employed 
in our analysis are polynomials constructed from com¬ 
binations of the unity, and s^. Therefore, in our 
FESRs, the leading contributions from the OPE arise 
solely from Cg and Cg. We neglect subleading logarith¬ 
mic corrections and treat all as constants. 

The DV contribution to the sum rules can be writ¬ 
ten as 

r“ ds nv 

(6) 

Jso ■*0 

where p°^(s) is the DV part of the spectral function in a 
given channel 

-Imn[^7\i), (7) 

TT 

that, for s large enough, we parametrize with the Ansatz 
of Ref. [18, 19] 

Pv/a(«) = exp(-(5v/A - Jvias) sin(av/A + Pvias) ■ (8) 

This adds 4 new parameters in each channel. 

We do not restrict our analysis to pinched weight 
functions. Rather, we include a weight function that is 
not pinched in order to constrain the DV parameters bet¬ 
ter. We work with three weight functions, namely, 

Wi = 1 , 

W2 = 1 - X^, 

W 3 = (1 - x2)(l +2x), (9) 

where x = s/sq. (The weight function W 3 is the one 
determined by kinematics and that yields Rj.) This 

’ We have checked that explicitly. 

^An alternative view on the dimension 2 contribution can be found 
in Refs. [23, 24] 


choice is motivated by the fact that we want to perform a 
self-consistent analysis, including all leading order con- 
tribuitons in the OPE — without truncating this series 
arbitrarily. The extensive explorations of Refs. [7, 8 ] 
have shown that this set of weight functions fulfills these 
requirements and allows for a good determination of a^. 
We remark that these weight functions also have good 
perturbative behaviour, in the sense of the analysis per¬ 
formed in [15]. 

3. Fits 

We have performed several different fits, combining 
subsets of the moments of Eq. (9), fitting to the V or 
the combined V and A channels. In the fits we include 
many different values of sq that lie inside a window 
[■Smin, Smax] in which our treatment of the perturbative 
series, the OPE contributions, and the DVs give an accu¬ 
rate description of the QCD correlator. One has to vary 
the value of Smin in the fits to check for stability. Fits 
to moments of a single weight function are performed 
minimizing a standard Fits that involve moments of 
more than one weight function, on the other hand, have 
too strong correlations to allow for a fit of this type. In 
this case, one must resort to other measures of fit quality 
and change the error propagation accordingly in order 
to account for the strong correlations. The procedure 
we adopt is discussed in detail in Appendix of Ref. [7]. 

We performed a number of consistency tests to assess 
the robustness of the outcome of our fits. To corroborate 
the results, we performed a study of the posterior prob¬ 
ability with a Markov-Chain Monte Carlo. The results 
were also checked for consistency comparing fits to V 
and combined analyses of V and A channels. Another 
important test is the stability against variations in sq, 
the upper limit of integration in Eq. (2). The descrip¬ 
tion of moments of the spectral functions must be valid 
for s < m^. We have checked this stability for several 
weight functions. Finally, the outcome of our fits was 
checked for consistency using the Weinberg sum rules 
(see Sec. 4.2). 

4. Main results 

4.1. Results for Os 

The results for that we obtain from the different 
fit set-ups are consistent within statistical errors. We 
choose to quote as our final value the one obtained from 
a fit to the V channel combining moments of the three 
weight functions wi, W 2 , and W 3 , of Eq. (9). Fits includ¬ 
ing the A channel require an extra assumption, namely. 
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that the asymptotic regime assumed in Eq. (8) has al¬ 
ready been reached although one works on the tail of the 
a I resonance. Although the consistency of the results 
indicate that this assumption may be fulfilled we pre¬ 
fer to rely on results from V channel only when quoting 
our final as values. Our final asinil) values for Nf = 3 
within the MS scheme are 

asiml) = 0-296 + 0.010 (ALEPH, EOPT), 
as(ml) = 0.310 + 0.014 (ALEPH, CIPT). 

The errors are dominated by statistics, but they also in¬ 
clude an estimate of the error due to varying the so win¬ 
dow and the error due to truncation of the perturbative 
series. When evolved to these results read 

as(ml) = 0.1155 + 0.0014 (ALEPH, EOPT), 

as(ml) = 0.1174 + 0.0019 (ALEPH, CIPT). 

These values are compatible with the ones obtained 
from the OPAL data [8]. Since the data sets are inde¬ 
pendent, the as values are virtually uncorrelated. This 
allows for a weighted average to be performed. The 
averaged values are slightly higher, since values from 
OPAL data tend to be larger. We find 

a,(m?) = 0.303 + 0.009 (ALEPH and OPAL, EOPT), 
asiml) = 0.319 + 0.012 (ALEPH and OPAL, CIPT), 

and at the Z boson mass one has 

as(ml) = 0.1165 + 0.0012 (ALEPH and OPAL, EOPT), 
as(ml) = 0.1185 + 0.0015 (ALEPH and OPAL, CIPT). 

4.2. Weinberg sum rules 

Results of combined fits to sum rules of the V and A 
channels allow for tests of the Weinberg sum rules. The 
first and second Weinberg sum rules can be written as 

X CXD 

ds(py-pl>)-2f^^0, (10) 

X oo 

dss{py-dy)-2mlf^^0, (11) 

where the pion pole contribution has been separated. 
Lor the second sum rule we assumed that terms of or¬ 
der m, mj, with i, j = u, d, can be neglected. To check 
such sum rules, one defines a point Ssw below which one 
uses the ALEPH data to compute the integral, and above 
which one extrapolates the DV Ansatz of Eq. (8) with 
parameters obtained from a fit to V and A. To be con¬ 
crete, we take the parameter values of a fit to V and A 


sum rules constructed from the three weight functions 
of Eq. (9) within CIPT, shown in Tab. V of Ref. [22]. 
(The results are very similar if another fit set-up is cho¬ 
sen.) We find, for the sum rules of Eq. (10) and (11), the 
results shown in Eig. 1 as a function of the point Ssw 




Figure 1: Weinberg sum rules, Eqs. (10) and (11). Data is used for 
5 < Ssw and the DV Ansatz of Eq. (8) is used for s > i'swasing results 
from a fit to y and A with moments of wi, W 2 , and ws, within CIPT. 
The parameters are given in Tab. V of Ref. [22]. 

To illustrate the importance of the DVs we display 
in Eig. 2 the results for the Weinberg sum rules as a 
function of Ssw but without the DV contribution. The 
comparison between these two results gives us confi¬ 
dence that the DVs obtained from the data are sound 
since their extrapolation to infinity is in excellent agree¬ 
ment with the constraints from the Weinberg sum rules. 

5. Conclusions 

We have extracted and non-perturbative parame¬ 
ters from the updated ALEPH data for hadronic t de¬ 
cays. This extraction is sound and our results fulfill 
consistency tests. The as values can be averaged with 
those obtained from the OPAL data. The final results 
are given in Sec. 4.1. The ALEPH data constrain the 
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(a) First Weinberg sum nale, Eq. (10), without DVs. 



(b) Second Weinberg sum rule, Eq. (11), without DVs. 

Figure 2: Weinberg sum rules, Eqs. (10) and (11) without the DV 
contribuiton. These results should be compared with Fig. 1 where the 
DV contribution is taken into account. 

DV parameters better than OPAL’s and our DV param¬ 
eters fulfill the constraints imposed by Weinberg sum 
rules. Nevertheless, the non-perturbative physics lim¬ 
its the accuracy with which the strong coupling can be 
extracted. 
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